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Abstract.  The  accuracy  of  astrometric  observations  conducted  via  a 
space-borne  optical  interferometer  orbiting  the  Earth  is  expected  to  ap¬ 
proach  a  few  microarcseconds.  Data  processing  of  such  extremely  high- 
precision  measurements  requires  access  to  a  rigorous  relativistic  model  of 
light  ray  propagation  developed  in  the  framework  of  General  Relativity. 
The  data-processing  of  the  space  interferometric  observations  must  rely 
upon  the  theory  of  general-relativistic  transformations  between  the  space¬ 
craft,  geocentric,  and  solar  bary centric  reference  systems  allowing  unique 
and  unambiguous  interpretation  of  the  stellar  aberration  and  parallax 
effects.  On  the  other  hand,  the  algorithm  must  also  include  physically 
adequate  treatment  of  the  relativistic  effect  of  light  deflection  caused  by 
the  spherically-symmetric  (monopole-dependent)  part  of  the  gravitational 
held  of  the  Sun  and  planets  as  well  as  the  quadrupole-  and  spin-dependent 
counterparts  of  it.  In  some  particular  cases  the  gravitomagnetic  held  in¬ 
duced  by  the  translational  motion  of  the  Sun  and  planets  should  be  also 
taken  into  account  for  unambigious  prediction  of  the  light-ray  dehection 
angle.  In  the  present  paper  we  describe  the  corresponding  software  pro¬ 
gram  to  take  into  account  all  classical  (proper  motion,  parallax,  etc.)  and 
relativistic  (aberration,  dehection  of  light)  effects  up  to  the  microarcsec¬ 
ond  threshold  and  demonstrate,  using  numerical  simulations,  how  ob¬ 
servations  of  stars  and/or  quasars  conducted  on  board  a  space  optical 
interferometer  orbiting  the  Earth  can  be  processed  and  disentangled.  For 
numerical  simulations  the  spacecraft  orbital  parameters  and  the  telescope 
optical-system-characteristics  have  been  taken  to  be  similar  to  those  in 
the  Hipparcos  mission.  The  performed  numerical  data  analysis  verihes 
that  the  relativistic  algorithm  chosen  for  data  processing  is  convergent 
and  can  be  used  in  practice  to  determine  astronomical  coordinates  and 
proper  motions  of  stars  (quasars)  with  the  required  microarcsecond  pre¬ 
cision.  Results  shown  in  the  paper  have  been  obtained  with  a  rather 
small  number  of  stars  (a  few  thousand).  Simulations  based  on  a  much 
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larger  number  of  stars,  e.g.,  from  the  Guide  Star  Catalogue  used  to  model 
original  observations  will  give  more  complete  information  about  potential 
abilities  of  the  space  astrometric  missions. 


1.  Introduction 

Current  technological  achievements  have  greatly  improved  the  accuracy  of  op¬ 
tical  observations.  Nowadays  an  optical  interferometer  based  on  board  a  space¬ 
craft  orbiting  the  Earth  is  considered  to  be  the  most  promising  tool  for  fun¬ 
damental  astrometry  and  astrophysics.  Microarcsecond  (and,  probably,  sub- 
microarcsecond)  accuracy  in  measuring  angular  distances  between  pairs  of  point¬ 
like  celestial  objects  (stars,  quasars,  etc.)  in  the  sky  could  become  feasible.  In 
practice,  the  accuracy  of  the  planned  space  astrometric  missions  is  anticipated 
to  reach  50  gas  (Johnston,  2000)  through  4  gas  (Danner  and  Unwin,  1999) 
to  1  gas  (Mignard,  2000).  It  is  commonly  recognized  (Kovalevsky,  2000;  S  off  el, 
2000)  that  data  processing  of  optical  observations  attaining  that  accuracy  cannot 
be  tackled  in  the  framework  of  the  classical  approach  based  on  the  Newtonian 
concepts,  which  must  be  abandoned  and  replaced  with  one  based  on  general 
relativity.  Indeed,  the  post-post-Newtonian  deflection  of  light  caused  by  the 
Sun  as  well  as  post-Newtonian  deflections  of  the  light  ray  associated  with  the 
oblateness  of  the  Sun  and  major  planets,  and  their  translational  and  rotational 
motions  can  achieve  values  of  the  order  of  1-10  gas  (Brumberg  et  crl. ,  1990). 

The  fact  that  observations  are  carried  out  on  board  the  Earth-orbiting  satel¬ 
lite  requires  precise  definition  of  the  satellite,  geocentric,  and  solar  barycentric 
reference  frames  along  with  the  corresponding  relativistic  transformations  be¬ 
tween  them.  Combining  relativistic  description  of  the  light  ray  trajectory,  writ¬ 
ten  in  the  solar  barycentric  reference  frame,  with  the  relativistic  transforma¬ 
tions  between  the  different  frames  allows  us  to  develop  unambiguous  models 
of  stellar  aberration  (Klioner  and  Kopeikin,  1992),  annual  and  diurnal  paral¬ 
laxes  (Kopeikin  and  Schafer,  1999;  Klioner,  2000)  of  stars.  Unfortunately,  the 
complexity  of  this  mathematical  model  makes  it  difficult  to  estimate  the  real 
accuracy  of  measurements  of  classical  and  relativistic  effects.  Furthermore,  it  is 
not  self-evident  that  the  model  in  question  will  make  it  possible  to  disentangle 
billions  of  measurements  of  arcs  between  pairs  of  different  stars  in  order  to  get 
their  astronomical  positions,  proper  motions,  and  parallaxes.  For  these  reasons, 
one  needs  to  study  these  and  related  questions  using  Monte-Carlo  numerical  sim¬ 
ulations  of  the  space  interferometric  optical  data-processing  before  the  mission 
is  launched. 

For  the  numerical  study  of  the  problem  under  consideration  the  software 
called  Ephemerides  Research  in  Astronomy  (ERA)  and  developed  in  the  Insti¬ 
tute  of  Applied  Astronomy  (Krasinsky  and  Vasilyev,  1996)  has  been  applied. 
This  software  was  designed  to  provide  a  universal  method  of  data  processing  of 
different  types  of  astronomical  observations  —  radio,  optical,  X-ray,  etc.  —  in 
ephemeris  astronomy.  Unfortunately,  the  classical  model  of  the  data  reduction 
of  optical  astrometric  observations  adopted  in  ERA  did  not  take  into  account 
all  relativistic  effects  properly.  Therefore,  the  goal  of  the  present  work  is  to 
re-design  the  model  of  the  data  reduction  procedure  used  in  ERA  on  the  basis 
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of  a  more  advanced  relativistic  model.  In  this  paper  we,  first,  describe  the  main 
modifications  used  in  ERA  extending  its  applicability  to  reach  the  microarcsec- 
ond  threshold  and,  second,  apply  the  modified  ERA  for  a  numerical  study  of 
processing  optical  data  obtained  with  microarcsecond  accuracy. 


2.  Relativistic  model  for  data  processing  in  microarcsecond  space 
astrometry 

In  order  to  avoid  misinterpretation  of  observed  astrometric  data  made  in  space 
at  the  microarcsecond-level  of  accuracy,  a  general  relativistic  description  of  light 
propagation  from  the  source  of  light  to  observer  and  that  of  reference  frames 
in  the  solar  system  are  required.  The  most  complete  and  successful  theoretical 
study  of  this  problem  was  undertaken  in  a  series  of  our  publications  (Brum- 
berg  et  al. ,  1990;  Klioner  and  Kopeikin,  1992;  Kopeikin  and  Schafer,  1999; 
Kopeikin,  1997)  triggered  originally  by  the  launch  of  the  Hipparcos  (Requieme, 
1985)  space-mission  and  the  work  on  the  POINTS  (Chandler  and  Reasenberg, 
1990)  project. 

The  approach  under  discussion  is  based  on  the  post-Newtonian  theory  of 
astronomical  reference  systems  and  the  mathematical  technique  of  matched 
asymptotic  expansions  (Kopeikin,  1988).  Three  basic  coordinate  charts  are  con¬ 
structed  by  solving  the  Einstein  equations  iteratively  in  the  slow-motion  and 
weak-field  approximation  —  the  barycentric  reference  system  (BRS)  of  the  so¬ 
lar  system,  the  geocentric  reference  system  (GRS),  and  the  satellite  reference 
system  (SRS).  The  BRS  covers  the  whole  space-time  and  is  used  to  describe 
light  propagation  from  the  observed  celestial  object  to  the  observer  as  well  as 
to  study  the  motion  of  bodies  around  the  Sun.  The  GRS  covers  a  space  domain 
in  the  vicinity  of  the  world  line  of  the  Earth’s  center  of  mass  and  is  used  to 
trace  the  motion  of  the  Earth-orbiting  satellite  carrying  the  space  optical  obser¬ 
vatory.  The  SRS  represents  a  local  coordinate  system  of  the  observer  (optical 
interferometer)  placed  at  its  origin. 

The  algorithm  for  reduction  of  space  astrometric  observations  performed 
on  board  Earth-orbiting  satellite  would  involve  the  relativistic  transformation 
of  star  positions  observed  in  the  SRS  to  the  barycentric  unit  vector  of  the  star 
defined  by  the  extrapolation  of  the  light-ray  trajectory  from  the  star  to  the 
observer  to  the  past-null  infinity  of  the  asymptotically-flat,  Minkowskian  space- 
time  where  the  inertial  reference  frame  is  formally  defined1.  More  specifically, 
the  data  reduction  procedure  includes  the  following  steps  (Klioner  and  Kopeikin, 
1992): 

1.  Components  of  the  observed  unit  vector  of  a  star  s'  defined  at  the  origin  of 
the  SRS  are  transformed  to  those  of  the  unit  vector  p  defined  in  the  BRS 
using  the  relativistic  transformations  from  the  SRS  to  the  GRS  and,  then, 
from  the  GRS  to  the  BRS.  This  two-step  transformation  eliminates  stellar 


JThe  advanced  post-Minkowskian  theory  of  light  propagation  including  explanation  of  the  rel¬ 
ativistic  terminology,  used  for  the  description  of  microarcsecond  astrometric  observations,  has 
been  recently  constructed  in  (Klioner,  2000)  (see  also  the  paper  by  Kopeikin  &  Gwinn  in  this 
proceedings). 
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aberration  caused  by  the  motion  of  the  satellite  about  the  Earth  and  the 
Earth  about  the  Sun. 

2.  The  positional  unit  vector  p  is  attached  to  the  point  of  observation  and 
still  includes  relativistic  deflection  of  light  due  to  the  gravitational  held  of 
the  solar  system.  It  is  eliminated  by  means  of  mapping  p  to  the  asym- 
potically  hat  space-time  using  the  solution  of  the  equations  of  light  prop¬ 
agation.  This  procedure  transforms  the  positional  vector  p  to  the  vector 
k  lying  in  the  hat  space-time  and  attached  to  the  point  of  observation. 
Microarcsecond  accuracy  of  observations  requires  us,  at  this  step,  to  ac¬ 
count  for  both  the  post-Newtonian  (1PN)  and  the  post-post-Newtonian 
(2PN)  perturbations  in  the  star’s  position.  More  specifically,  in  the  1PN 
approximation  dehections  of  hght  due  to  the  spherically-symmetric  part  of 
the  gravitational  held  of  the  Sun,  Moon,  Earth,  and  major  planets  must 
be  subtracted  along  with  the  dehections  caused  by  the  orbital  velocity  of 
Jupiter  and  quadrupole  component  of  the  held  of  the  major  planets  and 
the  Sun.  The  2PN  light  dehection  is  noticeable  only  for  the  gravitational 
held  of  the  Sun. 

3.  The  last  step  of  the  reduction  procedure  involves  taking  into  account  the 
parallactic  dispacement  and  proper  motion  of  the  star.  This  is  done  by 
making  an  ordinary  Taylor  expansion  of  the  vector  k  with  respect  to  the 
parallax  7r  and  time  t  —  to  increments  where  t  is  the  time  of  observation 
and  to  is  the  initial  epoch  of  observations.  In  this  way  we  obtain  the 
astrometric  position,  proper  motion,  and  parallax  of  the  star  which  can  be 
further  used  for  the  catalog’s  composition. 

Let  us  now  describe  how  the  model  works  in  the  analysis  of  numerically  simulated 

observations. 


3.  Numerical  simulation  of  astrometric  observations  and  their  pro¬ 
cessing 

The  modihed  program  package  ERA  (Krasinsky  and  Vasilyev,  1996)  was  used 
as  the  main  tool  for  the  numerical  simulation  of  observed  positions  of  celestial 
objects.  The  most  significant  modifications  have  been  done  in  that  part  of  the 
program  dealing  with  the  calculation  of  relativistic  hght-ray  deflections  and  the 
SRS-GRS-BRS  space-time  transformations.  The  new  program  has  been  care¬ 
fully  tested  by  calculating  residuals  in  star  positions  with  respect  to  the  previous 
version  of  ERA.  The  residuals  are  in  excellent  agreement  with  theoretically  pre¬ 
dicted  deviations  between  the  old  and  new  versions  of  ERA. 

To  simulate  the  space  astrometric  experiment  a  list  of  13000  stars  having 
random  proper  motions  and  radial  velocities  uniformly  distributed  over  the  ce¬ 
lestial  sphere  was  generated  and  treated  as  the  input  stellar  catalog  in  which 
positions,  proper  motions  and  radial  velocities  of  the  stars  were  considered  as 
precisely  known2.  It  was  also  assumed  that  all  stars  in  the  catalog  have  bright  - 


2We  did  not  include  parallaxes  of  stars  in  the  catalog.  They  are  assumed  to  be  taken  into 
account  in  the  future,  more  enhanced  study. 
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nesses  that  are  sufficient  to  reach  a  microarcsecond  observational  precision.  A 
fictitious  observer  was  placed  at  the  geostationary  satellite  carrying  three  tele¬ 
scopes  with  optical  axes  lying  in  one  plane  (the  satellite  “equatorial”  plane)  and 
separated  at  angles  54°,  78°  and  132°.  Each  telescope  was  assumed  to  have  a 
strip-like  coherent  held  of  observation  1.6°  wide.  The  satellite  is  assumed  to  be 
rotating  around  the  spin  axis,  perpendicular  to  the  “equatorial”  plane,  according 
to  a  scanning  law  that  is  similar  to  that  used  in  the  Hipparcos  mission3 

The  simulated  observations  generating  code  applies  the  procedure  described 
in  Section  2  in  the  reverse  order  to  simulate  positions  of  stars  as  they  would  be 
visible  from  the  satellite.  The  random  error  of  a  single  “observation”  was  taken 
to  be  1/ias.  Then,  the  program  starts  to  rotate  the  satellite  at  the  prescribed 
angular  velocity  and  “observations”  commence  at  some  epoch  to  such  that  the 
telescopes  sweep  over  the  simulated  positions  of  stars  step  by  step.  The  program 
calculates  the  telescope’s  directions  with  respect  to  SRS  and  looks  through  the 
input  catalog  reduced  to  the  SRS  frame  choosing  pairs  of  stars  which  are  avail¬ 
able  for  observations  simultaneously  by  two  of  three  telescopes.  Each  “observed” 
pair  of  stars  gives  one  conditional  equation  to  determine  parameters  of  the  stars 
and  magnitude  of  measurable  classical  and  relativistic  effects.  Due  to  the  limited 
computer  resources,  the  simulation  of  observations  have  been  done  at  the  span  of 
one  year  with  a  step  of  0.1  day.  This  gave  us  approximately  3000  observed  arcs 
between  pairs  of  stars  on  the  sky.  For  the  subsequent  parameter  estimation  we 
selected  only  those  stars  in  the  arcs  which  had  more  than  5  “observations”.  The 
number  of  such  stars  reached  1000.  Their  “observed”  positions  were  processed 
by  the  least-squares  method  in  accordance  with  the  model  of  section  2.  The 
post-fit  astrometric  positions  of  stars  were  compared  with  those  from  the  initial 
input  catalog.  The  next  section  presents  the  results  of  the  applied  procedure 
along  with  a  brief  discussion. 


4.  Results  and  discussion 

The  basic  results  of  our  numerical  study  of  the  problem  of  data  processing  of 
observations  done  with  microarcsecond  accuracy  on  board  an  artificial  satellite 
orbiting  the  Earth  are  as  follows: 

•  Even  a  small  number  of  observations  with  the  precision  of  a  single  mea¬ 
surement  approaching  1  pas  will  make  it  possible  to  improve  current  as¬ 
trometric  positions,  proper  motions,  and  parallaxes  of  stars4.  Relativistic 
effects  can  then  also  be  measured  with  improved  accuracy. 

•  There  can  exist  strong  correlations  between  classic  astrometric  parameters 
of  stars  and  those  characterizing  relativistic  effects.  The  correlations  get 
weaker  when  more  observations  are  included  and/or  distribution  of  stars 


3The  spin  axis  makes  a  constant  angle  of  43°  with  respect  to  the  direction  to  the  Sun  and 
precesses  around  this  direction  at  the  rate  of  6.4  rev/yr.  The  satellite’s  period  of  rotation 
around  the  spin  axis  was  assumed  to  be  128  minutes. 

4In  our  numerical  simulations  the  post-fit  accuracy,  e.g.,  astrometric  positions  of  stars  is  slightly 
better  than  40  g as. 
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all  over  the  sky  tends  to  be  more  uniform.  In  addition,  it  is  worth  not¬ 
ing  that  the  number  of  relativistic  parameters  depends  on  the  relativistic 
model  one  uses  for  data  processing.  More  straightforward  models  giving  a 
smaller  number  of  relativistic  parameters  will  be  more  preferable  for  space 
astrometry. 

•  The  number  of  observations  and  the  observed  distribution  of  stars  may 
crucially  restrict  the  feasibility  of  microarcsecond  space  astrometry.  More 
powerful  computer  codes  and/or  much  more  realistic  numerical  simulations 
including  bilhons  of  stars  are  required  to  answer  the  question  whether 
microarcsecond  threshold  in  measuring  classic  and  relativistic  effects  can 
be  reached  (see,  e.g.,  de  Felice  et  al.  (1998) 
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